A global analysis of inclusive diffract ive cross sections at HERA 
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We describe the most recent data on the diffractive structure functions from the HI and ZEUS 
Collaborations at HERA using four models. First, a Pomeron Structure Function (PSF) model, 
in which the Pomeron is considered as an object with parton distribution functions. Then, the 
Bartels Ellis Kowalski Wiisthoff (BEKW) approach is discussed, assuming the simplest perturba- 
tive description of the Pomeron using a two-gluon ladder. A third approach, the Bialas Peschanski 
(BP) model, based on the dipole formalism is then described. Finally, we discuss the Golec-Biernat- 
Wiisthoff (GBW) saturation model which takes into account saturation effects. The best description 
of all avaible measurements can be achieved with either the PSF based model or the BEKW ap- 
proach. In particular, the BEKW prediction allows to include the highest /3 measurements, which 
are dominated by higher twists effects and provide an efficient and compact parametrisation of the 
diffractive cross section. The two other models also give a good description of cross section measure- 
ments at small x with a small number of parameters. The comparison of all predictions allows us to 
identify interesting differences in the behaviour of the effective pomeron intercept and in the shape 
of the longitudinal component of the diffractive structure functions. In this last part, we present 
some features that can be discriminated by new experimental measurements, completing the HERA 
program. 

I. INTRODUCTION 

After many years, new measurements of the diffractive structure function have been made at the ep collider HERA, 
where 27.5 GeV electrons or positrons collide with 820 or 920 GeV protons [1-4]. The purpose of this paper is to 
analyse all these data under three different theoretical approaches, which are described in section II. We also study 
the compatibility of all experimental data, by fitting all sets together for each model. In this introduction, we remind 
briefly the main features of the models, used in this paper to understand the diffractive interactions. 

In a frame where the incident proton is very fast, the diffractive reaction can be seen as the deep inelastic scattering 
(DIS) of a virtual photon on the proton target, with a very fast proton in the final state. We can therefore expect to 
probe partons in a very specific way, under the exchange of a color singlet state. Starting with the pioneering theoretical 
work of Ref . [5] , the idea of a point-like structure of the Pomeron exchange opens the way to the determination of its 
parton (quark and gluon) distributions, where the Pomeron point-like structure can be treated in a similar way as (and 
compared to) the proton one. Indeed, leading twist contributions to the proton diffractive structure functions can be 
defined by factorisation properties [6] in much the same way as for the full proton structure functions themselves. As 
such, they should obey DGLAP evolution equations [7], and tehrefore allow for perturbative predictions of their Q 2 
evolution. On a phenomenological ground, these Diffractive Parton Distribution Functions (DPDFs) are the basis of 
MC simulations like RAPGAP [8] and they give a comparison basis with hard diffraction at the Tevatron [9], where 
the factorisation properties are not expected to be valid. Moreover, the study of DPDFs is also a challenge for the 
discussion of different approaches and models where other than leading twist contributions can be present in hard 
diffraction. Indeed, a strong presumptive evidence exists that higher twist effects may be quite important in diffractive 
processes contrary to non-diffractive ones which do not require (at least at not too small Q 2 ) such contributions. In 
fact, there are models which incorporate non-negligible contributions from higher twist components, especially for 
relatively small masses of the diffractive system. One of our goals in this paper is to take into account this peculiarity 
of diffractive processes. 

It is also useful to look at ep scattering in a frame where the virtual photon moves very fast (for instance in the 
proton rest frame, where the 7* has a momentum of up to about 50 TeV at HERA). The virtual photon can fluctuate 
into a quark-antiquark pair, forming a small color dipole. Because of its large Lorentz boost, this virtual pair has a 
lifetime much longer than a typical strong interaction time. Since the interaction between the pair and the proton is 
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mediated by the strong interaction, diffractive events are possible. An advantage of studying diffraction in ep collisions 
is that, for sufficiently large photon virtuality Q 2 , the typical transverse dimensions of the dipole are small compared 
to the size of a hadron. The interaction between the quark and the antiquark, as well as the interaction of the pair 
with the proton, can be treated perturbatively. With decreasing Q 2 the color dipole becomes larger, and at very low 
Q 2 these interactions become so strong that a description in terms of quarks and gluons is no longer justified, and 
the diffractive reactions become very similar to those in hadron-hadron scattering. 

The outline of the paper is as follows. In section II, we introduce the different models used to fit the HERA data 
that are described in section III. In section IV we start the analysis of these data with the PSF based approach, leading 
to the concept of parton distribution functions. We extract DPDFs under several fit conditions and many details are 
discussed in appendix A and B. Then, in sections V to VII, we follow with fits of HERA data using dipole models, 
for which we can consider saturation effects. In particular, the case of the two-gluon exchange model (BEKW) is 
shown to be quite efficient to reproduce the data over a large kinematic range including the higher- twist part. Results 
are discussed in section VIII. All the approaches covered in this paper give different properties concerning, e.g., the 
Pomeron intercept dependences and the longitudinal diffractive structure function. These aspects are discussed and 
and outlook is given. 

II. MODELS 

In this section, we discuss the phenomenology of the different models used to fit the HERA data. 

A. Pomeron Structure Function (PSF) model 

1. Reggeon and Pomeron contributions 

The diffractive structure function is investigated in the framework of Regge phenomenology and expressed 

as a sum of two factorized contributions corresponding to a Pomeron and secondary Reggeon trajectories [9]: 

if (3) (Q 2 ,/3,* P ) = /p /p (^p)F 2 i?(p) (Q 2 , / 3) + / R/p (*p)F 2 I?(R) (Q 2 ,/3) , (I) 

where xp is the fractional momentum loss of the incident proton and (3 has the form of a Bjorken variable defined 
with respect to the momentum P — P' lost by the initial proton as illustrated on Fig. 1 . 




FIG. 1. Schematic diagram of inclusive diffractive DIS, ep —* eXp. Four-momenta are indicated in parentheses. 



1 In a partonic interpretation, (3 is the momentum fraction of the struck quark with respect to the exchanged momentum 
P — P' (the allowed kinematical range of (3 is between and 1). 
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In the parametrisation defined by Eq. ( 1), F 2 can be interpreted as the Pomeron structure function and 

as an effective Reggeon structure function, with the restriction that it takes into account various secondary 
Rcgge contributions which can hardly be separated. The Pomeron and Reggeon fluxes are assumed to follow a Regge 
behaviour with linear trajectories av,n{t) = ap,]R.(0) + a P R t, such that 

fp/p,U/p(xp) = / 2a F . m (t)-l d ^ ( 2 ) 

Jt cut Xp 

where |t m »n| is the minimum kinematically allowed value of \t\ and t cu t = — 1 GeV 2 is the limit of the measurement. 
The values of the t-slopes parameters can be taken from Ref. [1] : a p = 0.06 GcV~ 2 , a B = 0.30 GcV~ 2 , Bp = 5.5 
GeV -2 and Br = 1.6 GeV -2 . Following the procedure described in Ref. [9], the Pomeron intercept ap(0) is let as a 
free parameter in the QCD fit and «r(0) is fixed to 0.50. 

The apparent Regge-factorisation breaking, in other words, the fact that we can not separate the ip from the 
(3 and Q 2 dependence of the diffractive structure function is explained in these kinds of models by introducing the 
Pomeron and Reggeon trajectories. As we see in the following, other models such as dipole models show an intrinsic 
factorisation breaking and it is of great interest to see if the secondary trajectories arc still needed by the fits. 



2. Extraction of partem densities in the Pomeron 

We assign parton distribution functions to the Pomeron and to the Reggeon. A simple prescription is adopted in 
which the parton distributions of both the Pomeron and the Reggeon are parametrised in terms of non-perturbative 
input distributions at some low scale Q 2 ,. For the structure of the sub- leading Reggeon trajectory, the pion structure 
function [11] is assumed with a free global normalization to be determined by the data. 

For the Pomeron, a quark flavour singlet distribution (zS(z, Q 2 ) = u + u + d + d + s + s) and a gluon distribution 
(zG(z,Q 2 )) are parametrised at an initial scale Qq, such that 

zS(z, Q 2 = Q 2 ) = [A s z Bs (1 - z) Cs (1 + D s z + E s Vz)] ■ , (3) 
zG(z,Q 2 = Q 2 )= [A G (l-zf°] -e^ , (4) 

where z — Xi/p is the fractional momentum of the Pomeron carried by the struck parton. In the following, we 
present the results in terms of quark density with the hypothesis that u = u = d = d= s = s, which means that any 
light quark density is equal to zS/6. 

Concerning the QCD fit technology, we fix as(M§) = 0.118 and the functions zS and zG of Eq. (3) and (4) arc 
evolved to higher Q 2 using the next-to- leading order DGLAP evolution equations with the code of Ref. [12]. All fits 
have also been performed with the code of Ref. [13] as a cross check which leads to the same results within 2%. The 
contribution to (/?, Q 2 ) from charm quarks is calculated in the fixed flavour scheme using the photon-gluon 

fusion prescription given in Ref. [14]. The contribution from heavier quarks is neglected. No momentum sum rule is 
imposed because of the theoretical uncertainty in specifying the normalization of the Pomeron or Reggeon fluxes and 
because it is not clear that such a sum rule is appropriate for the parton distributions of a virtual exchange. 

Note that diffractive distributions are process-independent functions. They appear not only in inclusive diffraction 
but also in other processes where diffractive hard-scattering factorisation holds. The cross section of such a process 
can be evaluated as the convolution of the relevant parton-level cross section with the DPDFs. For instance, the cross 
section for charm production in diffractive DIS can be calculated at leading order in a s from the j*g — > cc cross 
section and the diffractive gluon distribution. An analogous statement holds for jet production in diffractive DIS. 
Both processes have been analysed at next-to- leading order in a s . A natural question to ask is whether the DPDFs 
extracted at HERA can be used to describe hard diffractive processes like the production of jets, heavy quarks or 
weak gauge bosons in pp collisions at the Tevatron. The fraction of diffractive dijet events at CDF is a factor 3 to 
10 smaller than would be expected on the basis of the HERA data. The same type of discrepancy is consistently 
observed in all hard diffractive processes in pp events, see e.g. [15]. In general, while at HERA hard diffraction 
contributes a fraction of order 10% to the total cross section, it contributes only to about 1% at the Tevatron. This 
QCD-factorisation breaking observed in hadron-hadron scattering can be interpreted as a survival gap probability 
or a soft color interaction that need to be considered in such reactions, and a detailed computation of DPDFs from 
HERA can help in a better understanding of this phenomenon. 
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B. The Bartels Ellis Kowalski Wusthoff (BEKW) model 



This model assumes the simplest perturbative description of the Pomeron by a two-gluon ladder [16]. Contrary 
to the QCD fits described in the previous section, there is no concept of DPDFs in this approach. In Ref. [16], a 
parametrisation of the diffractive structure function in terms of three main contributions is proposed. The first term 
describes the diffractive production of a qq pair from a transversely polarised photon, the second one the production of 
a diffractive qqg system, and the third one the production of a qq component from a longitudinally polarised photon. 
We note that the higher Fock states such as qqgg are not included in this model as we discuss it further in the next 
section. More explicitely, we consider the modified form of the two-gluon exchange model, as expressed in Ref. [3], 
and labeled BEKW fits.' 

The three different contributions read 

ff»fe) = 4- /?(l-/3), (5) 

\XtpJ 

F 2 D(3) (qqg T ) = B as (ln[|j + 1]) (1 - (3) 1 , (6) 

= c (s) n4 {M + 7/4] ) V(i - wf ' (7) 

where 

«2(4) = "2(4) + "2(4) ln 

The ansatz for the /^-dependence is motivated by general features of QCD-parton model calculations: at small (3 
(large diffractive masses Mx) the spin 1/2 (quark) exchange in the qq production leads to a behavior ~ (3, whereas 
the spin 1 (gluon) exchange in the qqg term corresponds to (3°. For large (3 (small diffractive masses) pertubative 
QCD leads to 1 — (3 and (1 — f3)° for the transverse and longitudinal qq terms respectively. For the qqg term in qqg the 
situation is slightly more complicated : the exponent 7 is left as a free parameter. Concerning the Q 2 dependence, 
the first two terms are leading twist whereas the last one belongs to higher twist. The lnQ 2 -terms follow from QCD- 
calculations, and they indicate the beginning of the QCD Q 2 -evolution. Finally, the dependence on ip cannot be 
obtained from perturbative QCD and therefore is left free. An additional sub-leading trajectory [11] is added to the 
model, as for the DGLAP based fit, using the pion structure function to describe the Reggeon exchange. 

The free parameters in the fit are the three normalisation factors A, B and C, the i factors which describe 

the xjp dependence of the cross section and are not predicted by QCD, and the 7 parameter which describes the (3 
dependence of the qqg term. In the above prametrisation, n\ is set to zero, x is fixed to 0.1, Q 2 , is taken to be 1 
GeV 2 and as is taken at the fixed value of 0.25. The sub-leading trajectory [11] is added to the model, with an overall 
normalisation as a further free parameter, the flux definition being the same as in section II A 1. 



Qq 



(8) 



C. The Bialas Peschanski (BP) model 

The dipole model [19] provides another approach to describe the Pomeron as a perturbative object where high 
Fock states are introduced (qqgg, qqggg...). In this approach, two components contribute to the diffractive structure 
function [21]. First, an elastic component corresponds to the elastic interaction of two dipole configurations. It is 
expected to be dominant in the finite (3 region, for small relative masses of the diffractive system. Secondly, there is 
an inelastic component where the initial photon dipole configuration is diffractively dissociated in multi-dipole states 
by the target. This process is expected to be important at small (3 (large masses). 

A sub-leading trajectory [11] is added to the model, with an overall normalisation tir as a further free parameter, 
the flux definition being the same as in in section II A 1. The different contributions read 

pD{3) _ jyin pD(3),inel ^_ jyel pD(3) ,qel ^-el pD(3) ,qel ^ j^^pD{3) , Reggeon ^ 

The transverse and longitudinal clastic components can be written as 

pD(ei) = 12 N c e 2 a s 4 2ap+1 a3{xp) lQg3 _Q e _ o(xp) log 2 -3^ 
7r 2QoV/3 

(11) 



x/3(l-/3) 



^ [ -1^1 -p 



4 



p D(el) _ N c e 2 a s 4 2ap4 . 



.5, \i 2 Q -a(xp) log 2 n Q „ 

a (xp)log — — ;= e 2Q vTJ 



x/3 2 



n 2 



with 



and 



op = 1 + A 



l + ^41o g 2 



7T 1 

7C(3) a; 



The inelastic component is given by 



pD(in) _ 9 / z n 



1\ 7V c eV 



2Q - +1 a 3 ( 2 ;p)-^ e -^ log2 ^a5( / 3)/3- A 



(12) 
(13) 

(14) 

(15) 

(16) 



The free parameters used in the fit are ap(0), the exponents for the Pomeron and the secondary reggeon, the 
normalisations in front of the inelastic N m , transverse elastic iV|/ , longitudinal elastic Nf and the Reggeon component 
iVp, as well as the parameter Q , which is a typically non-perturbative proton scale. 

D. Golec-Biernat-Wiisthoff (GBW) saturation model 

The saturation model proposed some time ago by Golec-Biernat and Wiisthoff [22,23] was formulated in the color 
dipole picture. In this formalism both the inclusive and diffractive cross sections may be calculated. The diffractive 
structure function F^ 3>> is the sum of three contributions [23] : 



F 2 D(3) (Q 2 ,^,/3) = Ff + Ff + F™ S , 



(17) 



The incoming virtual photon 7* (transversely or longitudinally polarised) which interacts diffractively with the proton 
forms either the qq or qqg Fock state. The F^ 9 and F^ q components dominate in the region of large or intermediate 
values of (3 respectively, while the term F^? 9 is most important at small (3. The longitudinal component F™ is the 
higher twist contribution and therefore may be neglected. 

The two terms from Eq. (17) related to qq dipoles are given by the following formulae 



x P F° !L (Q 2 ,p,x P ) 



3Q 6 



32n 4 (3B D 



, /-1/2 

J^e 2 -2 / dzz 3 (l-z) 3 $ , 



(18) 



(19) 



where the sum runs over all active flavours /, each with charge e/ and mass m/, and integrate over the light-cone 
momentum fraction of the virtual photon carried by the quark or antiquark z. The lower limit of the integration 

is given by z = (1/2) ^1 — ^Jl — Am 2 , / Mj^j where Mx denotes the invariant mass of the diffractive system. The 

functions Q 2 and & i arc defined as 



$0,1 = 



Q 2 = z(l-z)Q 2 + m), 



/>oo 

/ rdrK ,i{Qr)a(x P ,r)Jo,i(kr) 
Jo 



(20) 
(21) 
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where a(xp,r) is the dipolc cross section specific for a given model. In the GBW model it has the form [22] 



a(x,r) = (7 {1 - cxp(-r 2 Q 2 at (x)/4)} 



) 2 

- sat 



w - (?)' 



(22) 



which introduces three parameters : the maximal possible value of the dipole cross section <tq and two parameters 
characterizing the saturation scale Q 2 at {x) that is A and x . 

The third component in Eq. (17) representing qqg state has the form 



x P F^JQ 2 ,p,x P ) 



81/3 as \—\ 9 



L 



(i-*)Q 2 



dkf In 



512n 5 B D 
(1 - z)Q'< 



dz 



z 



f 

Jo 



(1-Z)3 

udu a d (u/k t ,x XF )K 2 



1 - z 



(23) 



This formula was obtained in [23] using the two-gluon exchange approximation. In addition, strong ordering of 
transverse momenta of the quarks and the gluon in the qqg dipole was assumed. Notice that (23) is valid only when 
we consider massless quarks. 

Each contribution is proportional to the inverse of diffractive slope Bo for which we take Bp = 6.0 GeV -2 . The 
qqg component is also proportional to the strong coupling a s . Since this is not clear which value of a s should be 
used we made this coupling a parameter in our fits which effectively weights the qqg contribution to the diffractive 
structure function. 

In this paper we use the version of the GBW model with three light quarks only. In addition, we assume these 
quarks to be massless. 



III. DATA SETS 



In this section, we describe the data used in the following fits. We use the latest i-integrated diffractive cross 
section measurements from HI [1] (H1RAP) and ZEUS [3] (ZEUSMX) experiments, derived from the diffractive 
process e + p — > e + X + Y, where the proton stays either intact or in a low mass state Y. These two data sets are 
fitted independently and compared. 

In case of HI data [1] , measurements are presented in terms of the ^-integrated reduced cross section (xp, x, Q 2 ) 
defined by the relation 

d 3 a ep^eXY 4 2 2 

dxV^-^ {1 - y + Y^ { ix ^ Q) 

with (Jr^ = — i + ( y l _ y yz F®^\ We notice that the relation — F®^ is a very good approximation except 

at large y. 

ZEUS data are presented in terms of diffractive cross sections differential in Mx, da 1 di j f /dMx 1 and converted to 

values for the diffractive structure function values F®^ [3]. 

Note that these two data sets represent different methods of measurements, and therefore different domain in 
My, namely My < 1.6 GeV in case of HI and My < 2.3 GeV in case of ZEUS. In the following, all cross section 
measurements are corrected to the domain My < 1.6 GeV, after multiplying ZEUS values by the global factor 
0.85 [1,3]. 

In addition to the two previous data sets, we use the diffractive cross sections extracted by HI and ZEUS experiments 
from the process e + p — > e + X +p, where the leading proton is detected [2,4]. This way of tagging diffractive events 
is of course very interesting but suffers from a limited detection acceptance, which gives more restricted samples in 
terms of kinematic coverage. This is why we use them only in global fits with all available data sets. Also, in order 
to keep the compatibility between all data sets, cross section measurements from the leading proton samples are 
corrected to the domain My < 1.6 GeV, by multiplying them with the global factor 1.23 [1] which has been obtained 
by comparing the HI data requiring a rapidity gap and the data when the proton is tagged in the final state in roman 
pot detectors. 
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To summarize, we give the different data sets used in the fits performed in the following 2 . 

• er? measured by the HI collaboration using the rapidity gap method [1] labeled H1RAP, this is the default data 
set which is not corrected further; 

• F 2 D measured by the ZEUS collaboration using the U M X method" [3] labeled ZEUSMX, data multiplied by a 
global factor 0.85. Note that the ZEUS experiment gives directly the values of the structure functions , i.e. 
reduced cross sections corrected from the contribution of Fj?, and we use these values in the following; 

• er,P measured by the HI collaboration when a proton is detected in roman pot detectors [2] labeled H1TAG, 
data multiplied by the global factor 1.23; 

• F® measured by the ZEUS collaboration when a proton is detected in roman pot detectors [4] labeled ZEUSTAG, 
data multiplied by the global factor 1.23. 

IV. RESULTS USING THE PSF MODEL 

Before getting into the results of the QCD fits, let us mention a few technical points how the fits are performed. 
We do not fit the initial scale Qq as in Ref. [1], such that we get the best % 2 value and thus the best fit, with a 
fixed number of parameters. We have chosen to increase the number of parameters untill a fit result stable under a 
variation of the input scale Qq is obtained. As a result, the x 2 value and the shapes of the distributions for quarks 
and gluons do not depend on the choice of Qq. We have checked that we obtain very similar results as in Ref. [1] 
while doing the QCD fit under the same conditions as in Ref. [1] (see appendix A). 

As mentioned in the previous section, we consider first the HI data set alone (H1RAP). Further restrictions on the 
kinematical plane are requested so that the DGLAP based fits are valid. The two cuts, Mx > 2 GeV and (3 < 0.8, 
are necessary to avoid a kinematic range where higher twists effects are expected to be large, which would spoil the 
quality of the determination of the DPDFs. 

To allow the fit to be stable for the HI rapidity gap data alone (H1RAP), we have to consider 5 parameters for 
the quarks (see Eq. (3)) and 2 for the gluon (see Eq. (4)). With a similar cut scenario as in Ref. [1], Q 2 > 8.5 GeV 2 , 
M x > 2 GeV and (3 < 0.8, we obtain a X 2 = 169.3 for Ql = 3 GeV 2 and X 2 = 169.8 for Ql = 1.75 GeV 2 , for 190 data 
points included into the QCD fit. Only statistical and uncorrelated systematic uncertainties, added in quadrature, 
are considered in this analysis. We have also varied the values of Q 2 minl and the results are detailed in Appendix B. 
Following our method, we can produce DPDFs stable within variations of Q 2 min , which is not the case in Ref. [1]. It 
allows us also to fit all data points with Q 2 > 4.5 GeV 2 , whereas the cut is taken at 8.5 GeV 2 in Ref. [1]. In the 
kinematic domain defined by Q 2 > 4.5 GeV 2 , M x > 2 GeV and (3 < 0.8, we obtain a x 2 = 250.3 for 240 data points, 
using Qq = 3 GeV 2 as the initial scale for DPDFs definition. 

A parallel analysis is followed for the ZEUS data set [3] (ZEUSMX), with similar results illustrated on Fig. 15. Note 
that ZEUS measurements have been multiplied by the global factor 0.85 as explained in section III. The stability of 
the QCD fit is obtained with 4 parameters for the quarks (E$ = in Eq. (3)) and 2 for the gluon. In this case, x 2 
values and shapes are stable within variations of Qq and no significant dependence of the distributions is observed 
as a function of Q 2 min - In the kinematic domain defined by Q 2 > 4.5 GeV 2 , Mx > 2 GeV and < 0.8, we obtain 
a x 2 — 100.9 for 102 points and Qq = 3 GeV 2 , with the parameters given in Table I. For Q\ = 1.75 GeV 2 , we get 
X 2 = 99.9. Again for this study, only statistical errors are considered in this analysis. 



2 Note that not all models are able to describe the full data sets since their domain of validity is restricted and these additional 
cuts will be given while describing the fits 
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Parameters 


HI RAP 


ZEUSMX 


All data sets 


n' 2 


o vje v 


O vjc V 


O Vjrc V 


n 2 


A 5 PoV 2 
. J Vje v 


A * n P V 2 

'i.O Vjrtj V 


A Cn\l 2 


ajp 


i 1 on -t- n nn7 

1. 1ZU ± U.UU 1 


1.1U4 zt U.UUO 


1 lie ^ n nn^ 

1.118 ± U.UUO 


A a 








D 


ri iq i n no 
U.lo ± U.Uo 




n i /i -t- nil 

U.14 zt U.ll 


(la 


f) QO 1 A AO 


50 + Ofi 




D s 


6.14 ± 0.82 


5.65 ± 1.24 


5.72 ± 0.92 


E,s 


-3.98 ± 0.22 




-3.66 ± 0.19 


A G 


0.24 ± 0.06 


0.74 ± 0.15 


0.20 ± 0.06 




-fi 7fi + n iq 


Q Q« _|_ 1 1 « 


U. ( U —1— VJ . zi J- 


Ntr(HIRAP) 
N IR {ZEUSMX) 

N IR (H1TAG) 
N IR (ZEU STAG) 


5 77 ± 55 




6 65 ± 47 

5.06 ± 0.52 
4.64 ± 0.48 


X 2 1 (nb data points) 


250.3/240 


100.9/102 


377.6/444 



Table I- Pomeron quark and gluon densities parameters for the three different kinds of fits (H1RAP alone, ZEUSMX 
alone, and all four combined data sets). No Reggcon contribution is necessary for the ZEUSMX data set, as 
explained in Ref. [3]. Also, the parameter Bg is found to be for this particular fit. 



Data set 


X 2 (stat. 


and syst. 


error) 


Nb. of data points 


H1RAP 




217.9 




240 


HI TAG 




27.5 




57 


ZEUSMX 




109.5 




102 


ZEUSTAG 




22.7 




45 



Table II- x 2 values per data set for the global QCD fit with statistical and systematics errors added in quadrature 
(see text). The fit parameters are found to be very close if only statistical errors are used during the minimisation 

procedure. 

In a second step, we combine the four data sets defined in section III, with the normalisation factors for each set 
defined in this previous section. To realise this global fit, we use by default the total errors in the definition of the 
X 2 ■ The kinematic domain is again restricted to Q 2 > Q^im Mx > 2 GeV and (3 < 0.8, and we study the stability 
of the results when varying Qq and Q 2 min (see Appendix B). For Q 2 min = 4.5 GeV 2 , we get a x 2 value of 377.6 for 
444 data points, with the values for each data sets given in Table II and parameters listed in Table I 3 . Note that we 
assign a global Reggeon normalisation parameter to each data set for which this contribution is needed. However, the 
parameters for the Reggeon and Pomeron trajectories (defining the flux factors of Eq. (2)) are the same for all data 
sets, with ap(0) considered as a free parameter in the fit. 

The fit results are displayed in Fig. 2. We show the quark and gluon densities in the Pomeron (DPDFs) for HI 
data only (H1RAP), ZEUS data only (ZEUSMX) and all four data combined sets. We note that the Pomeron is 
gluon dominated for all fits. On the other hand, the shape of the gluon distribution is quite different between the HI 
(H1RAP) and ZEUS (ZEUSMX) data sets, where the high z contribution to the gluon density is small for the ZEUS 
data fit. The fact that the parton distributions for the combined fits are close to the HI results is simply due to the 
larger number of HI data points compared to ZEUS ones. The gluon density at high (3 is however poorly known and 
the uncertainty is of the order of 25 %. This high region is of particular interest for the LHC since it represents 
for instance a direct background to the search for exclusive events [9] . More data sets are needed to further constrain 
this region such as diffractive jets production at HERA, or Tevatron diffractive jet cross section measurements. 

It is also interesting to show the uncertainty on the high gluon density in the Pomeron by giving the uncertainty 
on the parameter v if the gluon density is multiplied by (1 — 0Y . The value of v from the fits to the H1RAP and 



3 We can also perform the QCD fits using the total error for each data set separately, which gives a x 2 value of 198.7 for 240 
data points in the case of H1RAP. When we compare to the global QCD analysis (with \ value of 377.6 for 444 data points), 
we observe that the \ 2 1 do] are compatible, which also justifies our approach to combine all data sets. 
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ZEUSMX data sets respectively is found to be ± 0.05 for a variation of % 2 of 15 units for the 4 combined data sets, 
reflecting the large uncertainty on the gluon density at large z. 



H1 data (H1 RAP) 

■■■■ ZEUS data (ZEUSMX) 
All data (4 combined sets) 




Z Z 



FIG. 2. Singlet and gluon distributions of the Pomeron as a function of z, the fractional momentum of the Pomeron carried 
by the struck parton, derived from QCD fits on H1RAP data alone, ZEUSMX data alone or the four data sets together. The 
parton densities are normalised to represent xp times the true parton densities multiplied by the flux factor at xp — 0.003. 



V. RESULTS USING THE BEKW MODEL 



Using formulae of section (II.B.l), we can perform fits of the structure function measurements (or reduced cross 
sections) over the whole kinematic range, including the highest /3 or lowest Mx values, as the longitudinal component 







of this model behaves as a higher-twist part. In a first step, a fit is performed on H1RAP data alone, with only one 
selection cut Q 2 > 4.5 GeV 2 : we obtain a \ 2 — 286.8 for for 247 data points. Only statistical and uncorrelated 
systematic uncertainties, added in quadrature, are considered. A parallel analysis is performed for the ZEUSMX data 
set, with a, x 2 ~ 193.5 for for 142 data points. Parameters are given in Table III for these two fits. 

Then, a global fit is performed to the reduced cross sections of all data sets, considering for this analysis the 
total errors (as in section IV). The parameters are given in Table IV and the values of the x 2 P er data sets arc 
presented in Table IV 4 . Fig. 4 shows the good description of the diffractive structure function measurements (H1RAP 
and ZEUSMX) by the fit result over the whole kinematic range. Only two data sets are presented on this Fig. 4, 

namely [1,3], to keep the plot lisible. From this Fig. 4, we observe that, taking into account the term F^ 3 \qq L ), the 
highest (3 bin is properly described except for the lowest Q 2 values. We notice also that the prediction for xpa^ ed 
is decreasing at low xp values for (3 — 0.9 : it is a reflection of the influence of the longitudinal component, which 
becomes large in this domain (see section III). If we apply the cuts Mj > 2 GeV and (3 < 0.8 to avoid this kinematic 
domain, the \ 2 f° r H1RAP improves considerably with a value of 226.6 for 240 data points and for H1TAG, we 
get a x 2 of 33.9 for 57 data points. For the other data sets (ZEUSMX and ZEUSTAG) the x 2 /dof is stable. The 
improvement for H1RAP comes mainly from the kinematic domain Q 2 < 15 GeV 2 and (3 = 0.9, which is not well 
described by the fit. 



Parameters 


HI RAP 


ZEUSMX 


All data sets 


A(qq T ) 
B{qqg T ) 
C(qq L ) 

n\ 

n\ 

n\ 

n\ 

7 

Xo 

Ql 


6.17 lO-^ ± 0.35 10-2 
0.72 10- 2 ± 0.04 10- 2 
9.55 10- 2 ± 5.94 10- 2 

0.09 ± 0.01 

0.04 ± 0.01 

0.35 ± 0.14 
0. 

9.45 ± 0.40 
0.1 
1. 


6.28 10- 2 ± 0.24 10-2 
0.70 10-2± 0.03 lO- 2 
7.59 10-2± 1.15 lO" 2 
0.18 ± 0.01 
0.01 ± 0.005 
0.10 ± 0.04 
0. 

9.32 ± 0.54 
0.1 
1. 


5.98 10-2 ± 0.26 lO" 2 
0.69 lO" 2 ± 0.03 10- 2 
6.57 10- 2 ± 1.30 10- 2 
0.12 ± 0.01 
0.04 ± 0.004 
0.15 ± 0.43 
0. 

9.76 ± 0.40 
0.1 
1. 


N IR (H1RAP) 
N IR (H1TAG) 
N IR {ZEU STAG) 


5.39 ± 0.57 




5.97 ± 0.49 
4.29 ± 0.53 
3.75 ± 0.50 


X 2 1 (nb data points) 


286.8/247 


193.5/142 


493.1/493 



Table III- Parameters for the BEKW fit performed in the range Q 2 > 4.5 GeV 2 . 



Data set 


X 2 (total error) 


Nb. of data points 


H1RAP 


308.3 


247 


HI TAG 


44.6 


59 


ZEUSMX 


116.5 


142 


ZEUSTAG 


23.7 


45 



Table IV- x 2 values (total errors) per data set for the global BEKW fit (see text). 

The resulting value of 7 is large (> 9), leading to a strong (3 dependence of the qqg term, which is dominant only in 
the low (3 region. In this kinematic domain, contributions in which more than two partons emerge from the diffractive 
hard scattering have also been found to be important in measurements of hadronic final state observables in diffractive 
DIS [17,18]. The longitudinal qq term becomes large for (3 > 0.5, indicating a strong higher twist contribution in this 
kinematic region. At intermediate (3, the transverse qq term, varying as (3(1 — (3) is dominant. These properties are 
illustrated in Fig. 3 for ip = 0.002. This shows that the scaling violations of diffractive structure functions observed 
experimentally can be described by the Q 2 dependence of the qq term at medium (3 since the exponent of xp depends 
on Q 2 . 



4 For this model, we can also perform the BEKW fits using the total error for each data set separately, which gives a \ 2 value 
of 226.2 for 247 data points in the case of H1RAP. When we compare to the global BEKW analysis (with \ 2 value of 493.1 for 
493 data points), we observe, as in the previous section, that the \ 2 /dof are compatible. 
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p p 




p p 



FIG. 3. The three different contributions of the two-ghion exchange model (BEKW), fitted on all data sets. Results are 
presented as a function of (3, for different Q 2 values, for a fixed ip = 0.002 value. 
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FIG. 4. Comparison of HI (H1RAP) and ZEUS (ZEUSMX) data sets with the prediction of the BEKW global fit. Only the 
statistical part of the uncertainty is shown for the data points on this plot. A dashed line is drawn for the prediction of the fit 
on points not included in the analysis. 
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VI. RESULTS USING THE BP MODEL 

In this section, we give the results of the fits based on the dipole model. Using formulae of section (II. B. 2), we can 
perform fits of the structure function measurements in the kinematic domain 4.5 < Q 2 < 120 GeV 2 , Mx > 2 GeV 
and j3 < 0.8. First, a fit is performed on H1RAP data alone, we obtain a \ 2 = 324.0 for for 175 data points. 
Only statistical and uncorrelated systematic uncertainties, added in quadrature, are considered. A parallel analysis 
is followed for the ZEUSMX data set, with a, \ 2 — 85.8 for 56 data points. Then, a global fit is performed to all data 
sets, considering the total errors. Parameters are given in Table V and the x 2 values per data set are shown in Table 
VI. Fig. 5 shows the result of the global fit compared to the data sets H1RAP and ZEUSMX. In the kinematic domain 
under study, the agreement is good but the extrapolation in the highest (5 — 0.9 fails to describe the data. We note 
also that in these fits, the Reggeon component is found to be very low, as the inelatic component is giving a similar 
behaviour. For the global fit, all Reggeon normalisations are found to be compatible with zero, but for H1RAP, for 
which wc get N R = 2.0 ± 0.5. 

The value of otp is found to be consistent with the expected intercept for a hard BFKL Pomeron [20]. This intercept 
is higher than the value obtained from the fit to the structure function Fi [19]. Qo is very close to the value obtained 
in the proton structure function fit. It should be noted that the scale Qo appears in a non-trivial way as the virtuality 
in the inelastic (Q/4Q ) and the elastic (Q/2v^9Qo) components. The validity of these results can be checked by 
starting with two different values of Qo for each component and the result of the fit leads to same values within errors, 
with the same \ 2 . Furthermore, imposing a constant scale for the elastic component, Q/2\//3Qo — * Q/Qo, leads to a 
very bad quality fit. 



Parameters 


HI RAP 


ZEUSMX 


All data sets 


ap 


1.407 ± 0.001 


1.329 ± 0.006 


1.339 ± 0.006 


Qo 


0.645 ± 0.005 


0.266 ± 0.010 


0.352 ± 0.012 


N m 


0.004 ± 0.0001 


0.003 ± 0.001 


0.003 ± 0.001 


N!jl 


28.105 ± 1.167 


154.11 ± 18.873 


124.230 ± 13.411 


N el 


1.129 ± 0.045 


33.904 ± 4.151 


27.330 ± 2.950 


X 2 1 (nb data points) 


324./175 


85.8/56 


601.3/478 



Table VII- Parameters for the BFKL fits (see text). 



Data set 


X 2 (total error) 


Nb. of data points 


H1RAP 


361.6 


232 


HI TAG 


32.8 


59 


ZEUSMX 


189.5 


142 


ZEUSTAG 


17.4 


45 



Table VIII- x 2 values (total errors) per data set for the global BFKL fit (see text). 
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FIG. 5. Comparison of HI (H1RAP) and ZEUS (ZEUSMX) data sets with the prediction of the BFKL global fit. Only the 
statistical part of the uncertainty is shown for the data points on this plot. A dashed line is drawn for the prediction of the fit 
on points not included in the analysis. 
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VII. RESULTS USING THE GBW SATURATION MODEL 



In this section, we give the results of the fits based on the dipole model including saturation effects. Using formulae 
of section (II. C), we can perform fits of the structure function measurements (or reduced cross sections) over the 
whole range of [3 and Mx values, as the longitudinal component of this model behaves as a higher-twist part. Since, 
as in other models, also in this case the sub-leading Reggeon part is added to the three terms in Eq. (17) the data for 
all values of xp are used in the fit. 

As shown in [24,25] the GBW saturation model works fairly well when fitted to the inclusive Fi data only when Q 2 
is smaller than 25 GeV 2 . This is because the Q 2 evolution which is expected to be important especially at high values 
of Q 2 is not present in the GBW model. This evolution introduced in [24] modifies the small r part of the dipole 
cross section to which the inclusive F 2 at high Q 2 is very sensitive. However the situation for diffractivc structure 
function F® is different since diffraction is sensitive rather to the intermediate r part of the dipole cross section which 
is essentially the same regardless of the presence or absence of Q 2 evolution. In fact, we have checked that imposing 
the Q 2 < 25 GeV 2 cut in the fit results in almost no improvement of the fit quality in terms of % 2 value. Hence, we 
decided to fit the data using only a Q 2 > 4.5 GeV 2 cut as for other models presented in this paper. 

First, a fit is performed on H1RAP data alone, we obtain a \ 2 = 272.9 for 247 data points. Only statistical and 
uncorrelated systematic uncertainties, added in quadrature, are considered in this analysis. A parallel analysis is 
followed for the ZEUSMX data set, with a x 2 ~ 268.0 for 142 data points. Then, a global fit is performed to all data 
sets and \ 2 = 564.5 is obtained for 493 data points. Parameters together with \ 2 values per number of data points 
are given in Table IX. Fig. 6 shows the result of the combined fit toghcter with the H1RAP and ZEUSMX data. A 
good description is found. 



Parameters 


H1RAP 


ZEUSMX 


All data sets 


<7 [mb] 


16.45 ± 0.24 


45.83 ± 0.25 


27.69 ± 0.14 


A 


0.239 ± 0.012 


0.16387 ± 


0.199 ± 


x 


6.78 • 10~ 3 ± 0.43 • 10~ 3 


1.59- 10- 8 ±0.02- 10~ 8 


2.539- 10" 5 ± 0.001 • 10~ 5 


a s 


0.131 ± 0.004 


0.075 ± 0.002 


0.097 ± 0.002 


N IR (H1RAP) 


4.60 ± 0.50 




5.02 ± 0.41 


N IR {ZEUSMX) 








N IR (H1TAG) 






2.96 ± 0.47 


Nm (ZEUSTAG) 






3.46 ± 0.51 


X^/(nb data points) 


272.9/247 


268.0/142 


564.5/493 



Table IX- Parameters for the GBW fits performed in the range Q 2 > 4.5 GeV 2 (see text). 



Let us recall that the original GBW parameters from [22] for the fit to F2 data with three light quarks are 
Co = 23.03 mb, A = 0.288 and xq = 3.04 • 10~ 4 . The value of a s used in [23] was fixed at a s = 0.2. Clearly the 
parameters obtained here for diffractive fits are substantially different from the original GBW parameters. The A 
parameter is smaller than 0.288 regardeless of the data set taken to the diffractive fit. The parameter xo varies a 
lot depending on which data we use. However, for the combined fit it is one order of magnitude smaller than found 
in [22]. This is also reflected in the form of the critical line, which is an important characteristic of the saturation 
model. This line in (x, Q 2 )-plane marks the transition to the domain where the saturation effect are important, i.e. 
to the left of it. In Fig. 7 we show the critical lines determined with the three sets of parameters from Table IX. Only 
the parameters which characterize the saturation scale, i.e. A and xq, are relevant in this case. The original result 
from [22] is also presented in Fig. 7 for reference. We see that the position of the critical line strongly depends on the 
data set used in the fit. The line obtained with the combined fit parameters is, however, not far away from the GBW 
result [22] whereas the line corresponding to ZEUSMX parameters is shifted significanly to the left with respect to 
other lines. 
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• ZEUS data (Mx) 
□ H1 data (LRG) 




FIG. 6. Comparison of HI (H1RAP) and ZEUS (ZEUSMX) data sets with the prediction of the GBW global fit. Only the 
statistical part of the uncertainty is shown for the data points on this plot. A dashed line is drawn for the prediction of the fit 
on points not included in the analysis. 
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Q 2 (GeV 2 ) 

FIG. 7. The critical line in the (x, Q 2 )-plane obtained with the three sets of parameters given in Table IX. The original GBW 
result from [23] is shown for reference. The shaded area shows the acceptance region of HERA. 
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VIII. DISCUSSION 



In the previous section, we have seen that all models give a correct description of all data sets. If we compare the 
X 2 /dof values per data sets in the kinematic range Q 2 > 4.5 GeV 2 , M x > 2 GeV and < 0.8, PSF and BEKW fits 
are compatible, with a global X 2 /dof of 0.83 and 0.81 (respectively). When we consider in addition the kinematic 
cut Q 2 < 120 GeV 2 , the BP and GBW approaches lead to a global x 2 /do,f of 1.24 and 1.31 (respectively). Then, 
PSF and BEKW predictions give a better description of the inclusive diffractive measurements but the BP and GBW 
models are also reasonnable and it makes sense to compare some specific features of all models considered above. 

Already, some properties have been mentionned in section (IV). For the BEKW parametrisation, we have seen that 
taking into account the term F^ 3 \qq L ) makes it possible to describe the highest bin (0 = 0.9), except for the 
lowest Q 2 values. Then, this model provides a very efficient description of all the inclusive diffractive measurements, 
with a simple functional form and few parameters. Of course, this is a leading order approach, as the dipole and 
saturation models we have considered in previous sections. The PSF model is the only one at ncxt-to-leading order, 
with the other advantage that DPDFs can be used to predict cross sections in hadron-hadron scattering, taking into 
account also the survival gap probability (see section II. A). 

Let us now study the dependence on xp which is directly related to the rapidity gap dependence. One can define 
an effective Pomcron intercept in the following way: 



where the t dependence is integrated out. This effective exponent can also be determinated for the inelastic and elastic 
components separately for the BP model. Results are presented for o?p as a function of xp (Fig. 8) or Q 2 (Fig. 9). 
In Fig. 8 (b), the effective intercepts of both components from the BP model and their sums have been presented. 
The range of obtained values sits essentially between these two limits except in the large xp region (xp > 10~ 2 ). It 
is clearly not consistent with the soft Pomeron value (1.08). The shape observed on Fig. 8 (b) can be explained by 
the large logarithmic corrections induced by the a 3 (xp) term, proportional to log 3 (l/xip), present in both diffractive 
components (see formulae 11, 13, 16). The effect of this logarithmic term induces also an xp dependence of the 
intercept. Moreover, it can be seen that the xp dependence of the intercept is different between the elastic and the 
inelastic components. This induces an apparent breaking of factorisation directly for the diffractive components of this 
model, which comes in addition to the known factorisation breaking due to secondary trajectories [21]. Fig. 9 shows 
the Q 2 dependence of the Pomeron intercept otp predicted by all models, with clear differences between the PSF fit 
and the BEKW model. In the case of PSF model, the Q 2 behaviour is by definition included in the DPDFs and not in 
the factorised xp part. In case of BEKW, there is a natural Q 2 dependence arising in the effective Pomeron intercept. 
A similar shape is observed in the case of the BFKL model, with still a normalisation factor compared to BEKW 
as discussed above. Note that the significant difference in behaviour between the PSF and the BEKW models is a 
challenge for experimentalists and the determination of the Q 2 dependence of the Pomeron intercept would provide 
a good discrimination between models. 

Another important conclusion from the previous analysis is the prediction for the longitudinal diffractive structure 
function by all models. Fig. 10 displays as a function of the longitudinal component of the proton diffractive 
structure functions (Fig. 10 (a)) and the ratio R of the longitudinal to the transverse components (Fig. 10 (b)). Here 
again, we observe significant differences. The BEKW of GBW models predicts only a longitudinal contribution at 
large 0, which behaves in Q 2 as a higher twist part (see section II). On the contrary, the PSF model predicts the 
leading twist part of F£ which is increasing at small values of and tends to when tends towards 1. It is 
interesting to observe that the BP model contains these two features for the behaviour of with a leading twist 
behaviour at low and a non-zero component for the largest values. 

Also, a striking feature of the diffractive proton structure functions concerns the scaling violation, namely the Q 2 
dependence at fixed xjp as a function of 0. The structure functions are increasing with Q 2 even at very high (see 
Fig. 11 and Fig. 12) at variance with the behaviour of the total proton structure function as a function of x. Note 
that we have only plotted the Pomeron part of the structure function, then subtracting the Reggeon component from 
the total structure function prediction. We observe that all models are in good agreement with some small differences 
in the highest bin, but much smaller than the experimental uncertainty accessible in this domain. 



We have discussed the most recent data on the diffractive structure functions from the HI and ZEUS Collaborations 
at HERA using four models, based on different approaches of the vacuum exchange structure. We have shown that 




(24) 



IX. CONCLUSION AND OUTLOOK 
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the best description of all avaible measurements can be achieved with either the PSF model or the BEKW model. In 
the case of the PSF approach, we have derived the DPDFs of the Pomeron under several variations in the analysis. In 
particular, a global analysis of all avaible data has been performed with a proper description of these measurements. 
The main features of the DPDFs, already observed in Ref. [9], are conserved. A large gluon content is still present 
with a large uncertainty of about 25% at large (3. Also, we have shown that the BEKW prediction allows to include 
the highest (3 measurements in the analysis. For these values of (3, higher twists effects are important and are correctly 
taken into account by the BEKW model through the longitudinal component of the diffractive cross section. We can 
mention that this model provides an efficient and compact parametrisation of the diffractive cross section over a large 
kinematic domain. The BP and GBW models also give a good description of all cross section measurements at small 
x with a small number of parameters. During the discussion section, we have noticed that the significant difference 
in behaviour for the effective Pomeron intercept between the PSF and the BEKW approaches is a challenge for 
experimentalists and the determination of the Q 2 dependence of the ap would provide a good discrimination between 
models. Concerning the longitudinal diffractive structure function, we have also observed interesting differences 
between models. In particular, the BP approach combines the low and high (3 behaviours of the PSF and BEKW 
models respectively. Once this observable could be accessed experimentally, it would be of great help for further 
developpements of the theory. 
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FIG. 9. Effective Pomeron intercept as a function of Q 2 . The dependence of ap on Q 2 is shown for all models at ip = 10 -3 
and (3 — 0.3 : in plot (a) we present predictions for BEKW as a full line, PSF fit as a dashed line, BP as a dotted line and 
GBW as a dashed-dotted line. A similar plot (b) is shown for the BP model alone, with elastic and inelastic component drawn 
separately. 
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FIG. 10. Predictions for x v F£ and R D = ^dKd as a function of f3 at Q 2 = 30 GeV 2 and x P = 10" 3 . We present 
predictions for BEKW as a full line, PSF fit as a dashed line, BP as a dotted line and GBW as a dashed-dotted line. 
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FIG. 11. Scaling violations, xpF 2 as a function of Q at different values of /3 and xp — 10 . We present predictions for 
BEKW as a full line, PSF fit as a dashed line, BP as a dotted line and GBW as a dashed-dotted line. 
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FIG. 12. Scaling violations, xpF 2 as a function of Q at different values of p and ip = 10 . We present the BP model 
predictions, with elastic and inelastic components shown separately (dashed and dotted curves respectively). 
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XI. APPENDIX A : REPRODUCING HI FIT A 



In Ref . [1] , the normalisations of the flux of Pomeron and Reggeon are defined in a different way, compared to our 
definitions in this paper or in our previous analysis (see Ref. [9] ) . In this appendix, to compare results from our fit to 
those of Ref. [1], we follow the same conventions as in Ref. [1] for the normalisations of the flux factors. Also, we use 
the Reggeon contribution from Ref. [11] and we perform a QCD fit under the same conditions as fit A of Ref. [1] : 
results are presented on Table A and are very close. Note also that, in our procedure, we use only the statistical and 
uncorrelated uncertainties added in quadrature, which gives a x 2 of 169.9 for 190 data points, whereas, in Ref. [1], 
the x 2 definition also includes parameters for all correlated systematics. 



parameters 


Our fit on HI data 


Table 3 of Ref. [1] (fit A) 


Q'o 


1.75 GeV 2 


1.75 GeV 2 


O 2 

^ mm 


8.5 GeV 2 


8.5 GeV 2 


«P 


1.118 ± 0.008 


1.118 ± 0.008 


As 


1.10 ± 0.32 


1.06 ± 0.32 


B s 


2.33 ± 0.35 


2.30 ± 0.36 


C S 


0.61 ± 0.15 


0.57 ± 0.15 


A G 


0.13 ± 0.03 


0.15 ± 0.03 


Cg 


-0.92 ± 0.16 


-0.95 ± 0.20 


N IR 


2.8 10~ a ± 0.4 10" 3 


1.7 10~ 3 ± 0.4 10" 3 



Table A- DPDFs parameters (see text). 
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XII. APPENDIX B : THE DGLAP FITS 

In this appendix, we want to give the systematics checks we performed on the DGLAP based QCD fits to the HI 
rapidity gap data alone (H1RAP), ZEUS M x data (ZEUSMX) or the combined data sets: 

• We check the dependence of the DPDFs on variations of the starting scale Qq in Fig. 13 for the HI data, very 
small changes are observed while changing the starting scale form 3 to 1.75 GeV 2 . Note that the results on the 
ZEUS Mx data and on the combined data sets are given in Figs. 15 and 16 and lead to the same conclusion 

• We check the fit stability by changing the cut on the lowest value of Q 2 of data to be included in the fit. 
The results are given in Figs. 14 (H1RAP), 15 (ZEUSMX) and 16 for the combined data set where we show the 
results of the fits to the HI and ZEUS data after applying a cut on Q 2 min of 4.5, 8.5 and 12 GeV 2 . Differences 
are noticeable at small (3 but well within the fit uncertainties. No systematic behaviour is observed within 
variations. 

• We check the dependence on the pion structure function by changing it from the Duke Owens parametrisation 
to the GRV one [10]. The results are given in Fig. 17 for the fits of all data sets and only very small differences 
are observed between both fits showing the weak dependence of the fit on the precise knowledge of the Reggcon 
structure function 

• We also check the fit stability by doing it with statistical errors only or with statistical and systematics errors 
added in quadrature. The results are shown in Fig. 18 and show very small different on the DPDFs between 
both fits. 
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FIG. 13. Singlet and gluon distributions of the Pomeron as a function of z, the fractional momentum of the Pomeron 
carried by the struck parton, derived from QCD fits on HI data. Results are presented with Q\ = 3 GeV 2 (full lines) and 
Qo = 1.75 GeV 2 (dashed lines). The parton densities are normalised to represent ip times the true parton densities multiplied 
by the flux factor at ip = 0.003. 
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FIG. 14. Singlet and gluon distributions ol the Pomeron as a function of z derived from QCD fits on HI data. Results are 
presented with Q 2 min = 4.5 GeV 2 (full lines), Q„ in = 8.5 GeV 2 (dashed lines) and Qmi„ = 12 GeV 2 (dotted lines). 
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ZEUS 

■■■■ ZEUS Q 2 min =8.5 GeV 2 
ZEUSQ 2 =1.75 GeV 2 




z 

FIG. 15. Singlet and gluon distributions of the Pomeron as a function of z derived from QCD fits on ZEUS data. Results 
are presented with Q^m = 4.5 GeV 2 (full lines), Q^in = 8.5 GeV 2 (dashed lines). Distributions obtained from a fit with 
Qq — 1.75 GeV 2 (dotted lines) is also shown. 
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All data 

■ ■ ■ ■ All data (Q 2 min = 8.5 GeV 2 ) 
All data (Q 2 = 1.75 GeV 2 ) 
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FIG. 16. Singlet and gluon distributions of the Pomeron as a function of z derived from QCD fits on all HI and ZEUS data 
sets. Results are presented with = 4.5 GeV 2 (full lines), = 8.5 GeV 2 (dashed lines). Distributions obtained from a 

fit with Qo = 1.75 GeV 2 (dotted lines) are also shown. 
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All data 

■■■■ All data (IR from GRV) 
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FIG. 17. Singlet and gluon distributions of the Pomeron as a function of z derived from QCD fits on all HI and ZEUS data 
sets. Results are presented for a QCD fit using the Owens (full lines) or GRV (dashed lines) parametrisations of the Reggeon 
contribution. 
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All data 

All data (stat. errors) 
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FIG. 18. Singlet and gluon distributions of the Pomeron as a function of z derived from QCD fits on all HI and ZEUS data 
sets. Results are presented for fits using the total errors (full lines) or the statistical errors (dashed lines). 
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